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We investigate gauge theories based on the smallest exceptional simple lie group G2. Our first 
model considered here is G2 Yang-Mills coupled to a fundamental Higgs field. In 4 spacetime 
dimensions we explore the phase diagram of the theory, showing that at larger Higgs masses 
the first order deconfinement phase transition turns into a crossover and therefore connects the 
low temperature confined phase with the high temperature deconfined phase. The second model 
investigated is G2 Yang-Mills coupled to fundamental fermions. It shares many features with 
QCD, especially fermionic baryons, but due to the absence of the fermion sign problem we can 
investigate this theory with Monte Carlo techniques even at low temperature and high baryonic 
density. First results on small lattices are presented. 
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1. Introduction 

Due to the fermion sign problem of QCD at real quark chemical potential, Monte Carlo techniques 
are not directly applicable to investigate the QCD phase diagram, especially at low temperatures 
and high densities [jl]]. Understanding the theory in this region of the QCD parameter space is very 
important for dense quark systems and the formation of compact stellar objects. Currently, the only 
reliable information is obtained from continuum methods, which usually require approximations 
or QCD-like theories having as many features as possible in common with QCD. One example of 
those theories is two colour two flavour QCD. Its phase diagram was explored in a series of papers 
[§, |3|]. The main drawback of two colour QCD is the absence of fermionic degrees of freedom in 
the hadron spectrum. In this work we investigate QCD-like theories based on the exceptional gauge 
group G2. It was already shown that G2 Yang Mills theory undergoes a first order phase transition 
from the low temperature confined phase to the high temperature deconfined phase [Q, ^. Addi- 
tionally, the gauge group SU (3) of strong interaction is a subgroup of G2 and this observation has 
interesting consequences [Q]. With a Higgs field in the fundamental 7-dimensional representation 
one can break the G2 gauge symmetry to the SU (3) symmetry via the Higgs mechanism. When the 
Higgs field in the action 

5[A,(/)] = J d^l^-^trF^^Ff"" + ^{Df,(j),D^(l)) + V{^)^ (1.1) 

picks up a vacuum expectation value v, the 8 gluons belonging to SU{3) remain massless and the 
additional 6 gauge bosons acquire a mass proportional to v. In the limit v — )• 00 they are removed 
from the spectrum such that G2 Yang-Mills-Higgs (YMH) theory reduces to SU (3) Yang-Mills 
theory. Even more interesting, for intermediate and large values of v, the G2 YMH-theory mimics 
SU (3) gauge theory with dynamical 'scalar quarks'. The second theory presented in this work is 
G2-QCD with the action given by 

S[A,y] = J d'x(^-^trFf,yF^''' + ^y{y^D^ +m + Yo^)vy (1-2) 

Here \j/isa massive Dirac spinor in the fundamental 7-dimensional representation of G2. We show 
that even at finite quark chemical potential /i the fermion determinant is non-negative. This allows 
us to investigate the phase diagram at zero temperature and finite chemical potential. Moreover 
the spectrum contains fermionic baryons and we expect the theory to behave in many aspects very 
similar to QCD. 

2. The confinement-deconfinement transition in G2 Yang-Mills 

For a gauge group with non-trivial center the trace of the Polyakov loop in representation ^ 

P^(x)=tr^^(x), =!^(x) = ^tr(^expi^\(T,x)t/T^, jSr = ^, (2.1) 

tranforms under center transformations like Pi^^Xx) i-> z*^P^(x) where z is an element of the center 
of the gauge group and k is the N-ality of the representation M. In a pure gauge theory the only 
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Figure 1: Left panel: Polyakov loop expectation values at the finite temperature confinement-deconfinement 
transition on a 16^ x 6 lattice. Right panel: Histograms of the Polyakov loop in the vicinity of the phase 
transition point Pciit — 9.765, pointing to a first order deconfinement transition. 



dynamical degrees of freedom are gluons in the adjoint center-blind representation. Therefore the 
Polyakov loop in any representation with non-zero N-ality is an order parameter for the spontaneous 
breaking of center symmetry. On the other hand the vacuum expectation value of the Polyakov 
loop is related to the free energy of an infinitely heavy test quark in the gluonic bath, (P) oc e^l^^ . 
Consequently, in the confined phase center symmetry is unbroken, while it is spontaneously broken 
in the deconfined phase. If we couple fundamental matter fields, as for example in QCD, the center 
symmetry is explicitly broken and the Polyakov loop is not an order parameter for confinement 
anymore. But it is still expected that QCD confines colour in the sense that charges of test quarks 
are screened by dynamical light quarks. In G2 gluodynamics the situation is very similar to the 
QCD case: Quarks and anti-quarks transform under the fundamental 7-dimensional representation. 
Since the center of G2 is trivial these fundamental charges can be screened by at least three gluons 
in the adjoint 14-dimensional representation, 

(7) (14) «) (14) (14) = (1) e . . . . (2.2) 

The Polyakov loop ceases to be an order parameter for confinement. In both, G2 gluodynamics 
and QCD, confinement can be identified with a linear rising inter-quark potential at intermediate 
scales and string breaking at larger distances in every representation of the gauge group [^, |8|]. 
Nevertheless on a finite lattice the Polyakov loop serves as an approximate order parameter which 
changes rapidly at the phase transition and is small (but non-zero) in the confining phase. We see 
a clear signal in the Polyakov loop at the confinement-deconfinement transition and the double 
peak structure in the histograms points to a first order confinement-deconfinement transition for G2 
gluodynamics (Fig. which is in agreement with earlier results in [0, ^. 

3. The G2 Higgs model 



The lattice action for the G2 Yang-Mills-Higgs theory reads 

SYMH[^,<J>]=i3£(l-^trRe^n) -jf£<I>,+^^,^<I>„ ^,-^, = 1, (3.1) 
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where is a seven component normalized real scalar field. In the limit jS ^ oo the gauge bosons, 
belonging to the coset space G2/SU (3) ~ S0{1) / S0{6) ~ S^, decouple and the theory reduces to 
an 5'0(7)-invariant nonlinear a-model. It shows spontaneous symmetry breaking down to S0{6) 
at a second order phase transition. With respect to the SU (3) subgroup of G2 the fundamental 
representations (7) and (14) branch into the following irreducible SU (3) -representations: 



(7)^(3)0(3)0(1) , (14)^(8)0(3)0(3). 



(3.2) 



S!'{7i) Gluodynamics 




Ci2 Gluodyni 



Figure 2: Expected phase diagram in the 
parameter space (7/^^, k). 



The Higgs field branches into a scalar quark, scalar anti- 
quark and singlet with respect to SU{3). Similarly, a 
Gi-gluon branches into a massless SU (3)-gluon and ad- 
ditional gauge bosons with respect to SU (3). The latter 
eat up the non-singlet scalar fields such that the spec- 
trum in the broken phase consists of 8 massless gluons, 
6 massive gauge bosons and one massive Higgs particle. 
In the limit of K" — )• 00 these 6 massive gauge bosons be- 
come infinitely heavy and the theory reduces to SU (3) 
Yang Mills theory. If K is lowered, in addition to the 8 
gluons of SU{3), the 6 additional gauge bosons of G2 
begin to participate in the dynamics. Similarly as dy- 
namical quarks and anit-quarks in QCD, they transform in the representations (3) and (3) of SU (3) 
and thus explicitly break the Z3 center symmetry. As in QCD they are expected to weaken the de- 
confinement phase transition [Q]. For k: = we recover G2 gluodynamics with a first order decon- 
finement phase transition. Fig. ^ shows the expected phase diagram. Our results for the complete 
phase diagram in the (j8, k) plane as calculated on 16^ x 6 up to 24^ x 6 lattices are summarized 
in Fig. |3[ We calculated histograms and susceptibilities of the Polyakov loop, the plaquette density 
and the Higgs action density near the marked points on the transition lines in this figure. From 
earlier calculations on smaller lattices we suggest that a triple point exists. An extrapolation to the 
point where the confined phase meets both deconfined phases leads to the couplings jStnp = 9.62(1) 
and K"trip = 1.455(5). Near this point the deconfinement transition is very weak, continuous or 
absent. Therefore we performed high-statistics simulations to investigate this region in parameter 
space more carefully. Up to a rather small region surrounding (jStdp, Kt^p) we can show that the 
deconfinement transition is first order and the Higgs transition is second order. But in this small 
window in parameter space, around the would-be triple point, the deconfinement transition is either 
second order or absent [|9|l. 



4. G2-QCD 



The lattice action for G2-QCD (1.2) reads 



SQCD[^,'i'] = j8 Ssy^m - -£»F.D^[^,m,At]»F- 



y ' 



(4.1) 



xy 



where Ssym is the Symanzik improved gauge action and D[^,m,/i] the Wilson-Dirac operator 
at bare fermion mass m and real quark chemical potential ^u. If a unitary operator T exists such 
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50(7) broken 
SU{3) deconfliiement 



50(7) unbroken 
G2 deconfinement 
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Figure 3: Phase transition lines on a 16 x 6 lattice. The solid line corresponds to the first order deconfine- 
ment transitions and the dashed line to the second order Higgs transitions. The plot on the right panel shows 
the details inside the small box in the plot on the left panel where the transition lines almost meet. The dotted 
line between the first order Unes corresponds to a window where the transition is a crossover 



that the Dirac operator obeys the relation D*T = T D with T* T = —1 and = T^^, the fermion 
determinant is non-negative [Q]. For Wilson fermions this implies the conditions T^T^ = ^* 
and Ty^T'^ = y^. Since every representation of G2 is real, we immediately find (for an Euclidean 
representation of the 7-matrices and charge conjugation matrix C) 

r = l^C75 ^ T*T = -1 detD['^,m,^]>0. (4.2) 

This feature of the fermion determinant makes Monte Carlo techniques applicable and allows us 
to investigate the complete phase diagram of G2-QCD on the lattice. The chiral symmetry of Nf 
flavor QCD is SU {Nf)L (E> SU {Nf)^ (S) U {1)b- Since the quark representation for G2 is real the 
dirac spinor decomposes into two Majorana spinors for vanishing mass. Consequently, we find 
an additional ^avowr symmetry and the continuous U{1)b turns into the discrete group Z(2). The 
chiral symmetry of G2-QCD is then [^] 

SUi2Nf)L=R^0Z{2)B. (4.3) 

Due to the discrete baryon number symmetry we can only distinguish between states with an even 
or odd number of quarks. Similar to QCD, mesons consist of an even number of quarks. Baryons 
are fermionic bound states with an odd number of quarks, for example, three quarks or a single 
quark and three gluons. Our first simulations are performed on rather small lattices with a spatial 
extend of L = 8 and T = 2 for finite temperature and on a 8"^ lattice in the case of zero temperature. 
The observables considered in this preliminary study are the Polyakov loop P, the chiral condensate 
X = and the quark number density Wq = y^^- First we compare the finite temperature 

phase transition at zero chemical potential to the phase transition in pure gauge theory. Our results 
for the Polyakov loop and the chiral condensate are shown in Fig. ^ (left panel). If we increase the 
hopping parameter from K' = Oto K" = 0.147, the transition in the Polyakov loop becomes weaker 
and the critical temperature decreases. This behaviour is also seen in the chiral condensate. The 
difference between its value in the chirally symmetric and broken phase grows with decreasing 
mass, due to the smaller explicit breaking of chiral symmetry. In Fig. ^ (right panel) we show the 
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Figure 4: Left panel: Polyakov loop (stars) and (rescaled) chiral condensate (dots) at finite temperature and 
zero chemical potential for ff = 0, 0.096, 0.131, 0.147. Right panel: Susceptibility of the Polyakov loop and 
chiral condensate. 




susceptibilites of both quantities. The peak becomes 
very broad for smaller quark masses, pointing to a cross- 
over. On these rather small lattices and within the given 
resolution in temperature, we can see no difference in 
the position of the peak. But whether the deconfinement 
and chiral transition critical temperature coincide and 
whether the transition really is a crossover at smaller 
quark masses has to be investigated in further simula- 
tions on larger lattices. At zero temperature the system 
stays in the vacuum until at approximately = 0.8 the 
quark number density increases. Fig. ^. At /i ?a 1.5 the 
system is saturated in the sence that each lattice site is 
occupied by an odd number of quarks. We also observe 
a transition in the Polyakov loop together with a drop in the chiral condensate. This points to a 
phase transition, or a crossover, from a confined phase to a chirally symmetric deconfined phase. 
At larger values of the chemical potential the Polyakov loop decreases again. This may be an lattice 
artifact due to the saturation of the system and has to be investigated further. If we increase the 
temperature again, the transition becomes weaker as shown in Fig. ^ 



Figure 5: Polyakov loop (stars), chiral con- 
densate (dots) and quark number density 
(triangles) at zero temperature. 



3.0 

2.5 - 
P 

ZO - 

1.5 - 

1.0 - 



•II, I j 



1.00 1.50 2.00 2.50 3.00 0.00 0.50 



[iiiilJtJJ 



loUMlli 

0.00 0.50 



1.50 2.00 2.50 3.00 



Figure 6: Polyakov loop P (left panel), chiral condensate x (center panel) and quark number density 
(right panel) for different temperatures. 
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5. Conclusions 

We have investigated QCD-like theories based on the gauge group G2. If we couple G2 Yang-Mills 
to a fundamental Higgs field, the theory mimics, in the broken phase, QCD with scalar quarks and 
anti-quarks with respect to the SU (3) subgroup of G2. In the limit of infinitely heavy quarks we 
recover SU (3) gauge theory with a first order deconfinement phase transition. If we decrease the 
mass the transition becomes weaker and, in a very small window in parameter space, turns into a 
crossover. This behaviour is already very similar to QCD. The second theory investigated is one 
flavour G2-QCD which has many features in common with QCD. In this first study we show that 
the finite temperature phase transition becomes weaker with decreasing mass of the quarks. At zero 
temperature and finite chemical potential we find a transition in the Polyakov loop and the chiral 
condensate together with an increasing baryon number density. First results on larger lattices will 
be published in a forthcoming paper [[l0|]. 
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